Zero sound in triplet-correlated superfluid neutron matter 

L. B. Leinson 
Abstract 

The linear response of a superfluid neutron liquid onto external vector field is studied for the 
case of 3 P2 — 3 -^2 pairing. In particular, we analyze the case of neutron condensation into the 
state with rrij = which is conventionally considered as the preferable one in the bulk matter of 
neutron stars. Consideration is limited to the case when the wave-length of a perturbation is large 
as compared to the coherence length in the superfluid matter and the transferred energy is small in 
comparison with the gap amplitude. The obtained results are used to analyse collisionless sound- 
like excitations of the superfluid condensate. Zero sound (if it exists) is found to be anisotropic 
and undergoes strong decrement below some temperature threshold depending substantially on the 
i ■ intensity of Fermi- liquid interactions. 
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I. INTRODUCTION 



At temperatures below the critical value T c the neutron Fermi liquid in the bulk matter 
of neutron stars is expected to develop a triplet superfluid condensate |l|, |2|]. Considerable 
work has done with the most realistic nuclearpotentials for determining the magnitude of the 



superfluid gap at different matter densities 



-|8j, while the low-energy collective excitations 



of such superfluid liquid are not well investigated as yet. In the meantime such excitations 
can play an important role in the evolution of neutron stars. For example, the decay of 
collective spin waves into neutrino pairs, occuring in a superfluid core of neutron star through 

neutral weak currents, presents a new mechanism of intensive cooling in some domain of 

S~[ I — I 

111] . A recent investigation [12] has shown sound-like oscillations of 

1 S superfluid neutron matter (called superfluid phonons) due to a very large mean- free- 
path influence heat conduction in a magnetized crust of neutron stars, where the motion 
of electrons is very anisotropic. Analogous effects might be expected in the neutron star 
core, since spontaneous breaking of rotation invariance and the baryon number caused by 
the triplet condensation should lead to the appearance of several Goldstone modes [131 ] . 

Generally speaking under /3-equilibrium the superdense core of neutron stars is composed 
of neutrons with a small admixture of protons and electrons. A fraction of hyperons can 
also appear at higher densities. It is well known however that long-range electromagnetic 
interactions push out collective oscillations of the charged particles up to the plasma fre- 
quency which is large sufficiently for (approximate) decoupling of the plasma modes from the 



sound-like oscillations of neutral component 14Ml6|. In this case the effect of short-range 



proton-neutron interactions is reduced mostly to renormalization of the effective mass of 
neutrons participating in the collective oscillations [171 . Il8| . Therefore the problem can be 
simplified considering the sound-like excitations in a pure neutron superfluid liquid. 

Previously the sound modes at finite temperatures have been investigated for isotropic 
singlet-spin superfluids [19] and for the case of triplet p-wave pairing in superfluid liquid 3 He 



2CH22J . Although a qualitative picture of the sound-like waves in superfluids is very similar, 
the above theories cannot be immediately applied to the case of superdense superfluid neu- 
tron matter. The well-developed theory of isotropic pairing cannot be applied because the 
triplet condensate in superdense nuclear matter is expected to be anisotropic. The theory of 
sound-like collective excitations in an anisotropic phase of superfluid 3 He is designed only 



for extremely fast waves with a velocity that is very large as compared to the Fermi velocity. 

Since spin-orbit and tensor interactions between neutrons are known to dominate at high 
densities the neutron pairing involves a mixing of 3 P2 and 3 F2 channels [3j-l8|. The sum 
of spin-orbit and tensor interactions cannot be described with the aid of a sole coupling 
constant. This complicates the investigation of collective excitations in standard ways with 
making use of an explicit form of the pairing interaction, where the sole coupling constant 
drops out of the equations by virtue of the gap equation. 

However, when the wave-length of the perturbation is large as compared to the coherence 
length in the superfluid matter and the transferred energy is small in comparison with the 
gap amplitude, as is typical for sound-like excitations, the collective motion of the condensate 
can be described in terms of total variable phase which (in the BCS approximation) can be 
derived immediately from the current conservation condition. This approach, for the first 



time suggested in Ref. 23], allows to avoid any explicit form of the interaction in the pairing 
channel. Residual Fermi-liquid interactions can be incorporated into the theory as a set of 
molecular fields [24| . In application to polarization functions, this approach is well developed 
in Ref. [13]. 

This paper is organized as follows. Section II contains some preliminary notes on how 
the Fermi-liquid interactions can be reduced to molecular fields. In Sec. Ill we derive, 
in the BCS approximation, the linear response of the triplet-correlated superfluid neutron 
liquid onto an effective vector field given by the sum of external and molecular fields. In 
Sect. IV we express self-consistently the effective fields via external fields, thus obtaining 
the linear medium response by taking into account the Fermi-liquid interactions. In Sect. V 
we analyze the poles of the longitudinal response function in order to derive the dispersion 
of sound-like oscillations in the condensate. Section VI contains a short summary of our 
findings and the conclusion. Throughout this paper, we use the system of units H = c = 1, 
and the Boltzmann constant ks — 1- 

II. FERMI-LIQUID INTERACTIONS AND MOLECULAR FIELDS 

It is well known that the Landau theory of a normal Fermi-liquid is based on the fact that 
a large part of the interactions can be taken into account with the aid of renormalizations 
effects. An effective Hamiltonian of the system contains the renormalized single-particle en- 



ergy of quasiparticles with occupation numbers n(p, er) and a residual interaction between 
changes 5n in the quasiparticle occupation at the Fermi surface. As has been shown by 



Leggett |25|], the quasiparticle pairing does not change the net occupation for a given direc- 
tion on the Fermi surface, if approximate particle-hole symmetry is maintained. Thus the 
Fermi-liquid interactions remain unchanged upon pairing. In other words, the Fermi-liquid 
interactions do not interfere with the pairing phenomenon. 

In our analysis we shall also assume that the anisotropy of the order parameter, which 
takes place in the case of triplet-spin pairing, plays no significant role in the Fermi-liquid 
interactions. This assumption is clearly justified because the characteristic length associated 
with the Fermi-liquid interactions is of the order of the inverse Fermi momentum, p^ 1 , and 
hence is much smaller than any other characteristic length entering the problem. This allows 
us to disregard the spin-dependent part of Fermi-liquid interactions in the vector channel 
we shall consider. 

Since we are interested in values of the neutron momenta near the Fermi surface, p ~ 
Pi?n, the amplitudes of the Fermi-liquid interactions / (nn') can be expanded into Legendre 
polynomials. In the vector channel these interactions are spin-independent and can be 
completely described in terms of the infinite set of Landau parameters Fi . In practice even 
for a saturated nuclear matter one does not know the Landau parameters Fi for I > 2, 
and in actual calculations they are frequently put equal zero. The remaining Fermi-liquid 
interactions can be written in the form 

gf(pri) = F + F 1 rm', (1) 

where g = ppM* /ir 2 is the density of states near the Fermi surface in the normal state, and 
the effective mass of a neutron quasiparticle is defined as M* = p F /vp, where vp is the 
Fermi velocity. 

This approach can be considered as a model of the Fermi-liquid interactions. It is known, 
however, that the Landau interactions with I > 2 do not modify the longitudinal response 
functions of normal (nonsuperfluid) Fermi liquid. For a one-component Fermi liquid this was 
demonstrated, for example, in Ref. 26| . The same result was obtained for a one-component 
1 S'o superfluid Fermi liquid in Ref. 27J where the effective vertices and the polarization 
functions (at u>, qv F <C A) have been found to depend only on F and Fi. 

Interactions (Q renormalize the normal energy of a quasiparticle in the weak external 



vector field A^ = (A , A) as 

e (p) ~ e (p) + c v A - ^pA + -F J2 Sn (p', a') + -\f iV ^ p'Sn (p', a) . (2) 

M e j^ QPf % 

We denote as e p = Vp (p — pp) the quasiparticle energy related to the Fermi energy in the 
normal state; cy is the coupling constant, which depends on the nature of the external field, 
and n (p', a') is the distribution function of neutron quasiparticles with momenta p and spin 
a. 

From Eq. (J2J) it is clearly seen that Fermi-liquid interactions can be reduced to molecular 
fields [2J], defined as 

A rnoi s _L Fq y- 6 n (p, a) , (3) 

CyQ 

and 

A^^-^L^pSnfav) . (4) 

c vB p F ^ 

Then Eq. (J2J) can be reduced to the form 

e(p)=e(p) + c v [A^-^A^) , (5) 

where the effective fields are given by the sum of external and molecular fields, 

M* 
A e ff = A + A™°\ A eff = (A + A mo/ ) . (6) 

The molecular fields depend on the charge perturbation and current density and should 
be calculated consistently. Therefore first we shall perform the calculation of the medium 
response onto the effective field and next find the explicit form of the effective fields in a 
self-consistent way. 

III. BCS RESPONSE IN THE LIMIT OF u, qv F «C A 

The triplet order parameter in the neutron superfluid is a symmetric matrix in spin space 
which can be written as 

V(n) = Ab(n)&g, (7) 

where a = (<7i, 0"2, 03) ar e Pauli spin matrices; g = i&2, with gg = — 1; and 1 is the 
2x2 unit matrix in spin space. In the ground state, the gap amplitude A is a constant 



(on the Fermi surface), and b (n) is a real vector in spin space which we normalize by 
the condition (b 2 (n)) = 1. Hereafter the angle brackets denote angle averages, (...) = 
(47r)~ J <in... . The angular dependence of the order parameter is represented by the unit 
vector n = p/p which defines the polar angles (8, <p) on the Fermi surface. In the components, 
ni = sin 9 cos cp, ri2 = sin 8 simp, n^ = cos 6. 

Making use of the adopted graphical notation for the ordinary and anomalous propa- 
gators, G = — ►— , G~(p) = — « — , F^ = » « , and F^ = « » , it is convenient to 
employ the Matsubara calculation technique for the system in thermal equilibrium. Then 
the analytic form of the propagators is as follows (see, e.g., Ref. |9[) 

G (r]n, p) = G (i] n , p) 5 a p , G~ (r] n , p) = G~ (r/ n , p) 5 aP , 

F« ( Vn , p) = F ( Vn , p) bag , F^ ( Vn , p) = F ( Vn , p) g&b , (8) 

where the scalar Green's functions are of the form G~ (r) n , p) = G (—i] n , — p) and 

G (^,P) = J" ^ , F(Vu,P) = 2 2 • (9) 

In ' p In ' p 

In the above, r] n = n (2n + 1) T with n — 0, ±1, ±2, ... is Matsubara's fermion frequency, and 
we assume the "unitary gap matrix" in the ground state, T> (n) T>* (n) oc 1, thus obtaining the 
energy of a one-particle excitation in the form E 2 = e 2 + A 2 b 2 (n), where the (temperature- 
dependent) energy gap, A n = A (T) b (n), is anisotropic. 

The following notation will be used below. We designate as L x> x (w, q; p) the analytical 
continuation of the Matsubara sums: 

Lxx< (w m ,p+-;p--J =T^X (rj n + uj m ,p+-)X' (r] n ,p--J , (10) 

n ^ 

where X, X' e G, F, G~, and u m = 2nTm with m = 0, ±1, ±2.... 

It is convenient to divide the integration over the momentum space into an integration 
over the solid angle dry and over the energy de p and operate with integrals 

Ixx' (w, n, q; T) = - / de p L xx , (u,p+-,p--J . (11) 

These are functions of u, q and the direction of the quasiparticle momentum p = pn. 
In deriving Eq. (TTT1) integration over de p is extended to — oo since the neutron matter is 
extremely degenerate. 



Consider the medium response onto the effective vector field ([6]) in the BCS approxima- 
tion. In this case the ordinary three-point vector vertices of a quasiparticle and a hole are 
defined in accordance with Eq. §5§: 

^ (p) = (1, p/M*) , 7; (P) = 1, (-P) • (12) 

We use greek letters for Dirac indices, /i = 0, 1, 2, 3. 

Variation of the anomalous self-energy T> (n) in the field A^f*, can be described with the 
aid of anomalous three-point vertices Tp, ' , defined as: 

c v fW = 8VM/8A£„ , c v fV = SD®/8AZ ff . (13) 

The anomalous vertices are 2x2 matrices in spin space which, near the Fermi surface, 
depend on the transferred energy momentum k^ = (u, q) and the direction n = p/p of the 
quasiparticle velocity 



The Ward identity implies the 
the order parameter in the system 



ollowing relations between the anomalous vertices and 



□ 



[28] (see also Refs. (l7|,[23|): 



k v T (i) = 2J) ( n ) ; PT^ 2) = -2£> f (n) . (14) 

We now restrict our consideration to the case, when the wave-length of the perturbation is 
large as compared to the coherence length and the transferred energy is small in comparison 
with the gap amplitude, u>, qvp ^ A. The only possible collective motion of the condensate 
in this case is a variation of the total phase without a change of the order parameter structure. 
Then the Ward identity reveals that for a uniform medium the anomalous vertices can be 
written in the form 

f« = Q, (u, q) bag , ff = -Q, (u, q) £<xb , (15) 

where the unknown vector function Q M (cu, q) satisfies the condition 

fc"Q„ = 2A . (16) 

For further progress, let us consider the retarded BCS polarization tensor in the vector 
channel Il^u^q). The latter can be found using the fact that the current in the system, 
j 11 = (8p,j) is connected to the linear correction 5G to the Green's function of a quasiparticle 

6 



in the effective external field A^** and can be obtained by analytic continuation of the 
following Matsubara sums 

5p = c v tY: f '0Tr (66) , (17) 

Pm 

J = TF-T Y I ^pTr (66) - ^A«" , (18) 

Pm 

where N = p F / (37r 2 ) is the total number density of neutrons. 

The linear correction to the Green's function of a quasiparticle caused by the external 
field A^ff is given by the diagrams shown in Fig. (TJ and can be written analytically as 



5G = » • » + »« • « * + 

FIG. 1. Correction to the ordinary propagator of a neutron quasiparticle in the external field. The 
ordinary and anomalous tree-point vertices are depicted by small circles and shadowed triangles, 
respectively. 

SG = c v (GG% - FF % (oh) (ah) 

+GF Tf>g (ah) + FG (ah) gf^) A% } , (19) 

where the anomalous vertices are to be taken in the form (fTo'j) . and we use the notation 
GG = G (r) n + u m , p + q/2) G (r) n , p - q/2), etc., 

Inserting Eq. (TT9l into Eqs. (ITTj) and ({TBI) we can derive the retarded polarization tensor 
II MV with the aid of the standard relation j M = IP^u;, c\)A e J^ . In this way after a little 
algebra we find 

n AM = Aa m + AaQ m , (20) 

where Q^ (u>, q) is still an unknown vector to be found. The functions A\^ (u>, q) and A\ (u>, q) 
are given by 

A AAt = c 2 v q I {i\i^T GG - jx%b 2 X FF ) + -^-6 X iS ilM J , (21) 

and 

A A = c vQ (2 lx b 2 l FG ) = c 2 vQ - < 7A (w + qv) b 2 X FF ) . (22) 



To obtain the second equality we made use of the identity 

1-fg — —%gf = — rr — Iff , (23) 

which can be verified by a straightforward calculation. 

The function Q^ (u, q) can be found from the requirement that the polarization tensor 
( |20l) must satisfy current conservation conditions, W" v k v = 0, and k u Il Ufl = 0, which can be 
written as the two coupled equations 

A^K + 2A A^ = , (24) 

fc„A v " + Kk v (q) Q" (?) = . (25) 

We made use of relation (TT6|) in the first equation. 

From the explicit form (J2"T|) for A^ u one can easily find the following relation 



g^A MA = A x ^ + c 2 v g ( ( 7a (w + qv) - (w - qv) 7a ) &%-*■) • (26) 

The first term on the right-hand side of this expression can be replaced as A^ v k v — > — 2A A M , 
in accordance with Eq. ( 1241) . Inserting the obtained result into Eq. ( 125|) . we arrive at 

2AA A - 4 e < ( 7A (w + qv) - (w - qv) 7a ~) 6 2 X ff ) 

Qa = ^v • ( 7) 

This formula can be simplified making use of the explicit form ( 1221) for A\. In the components 
we obtain the expressions: 

u (b 2 X FF ) , x 

<9o = 2A- ) ff^ , 28 

<(u; 2 -(qv) 2 )6 2 Z FF ) l J 

(v (qv) b 2 X FF ) , N 

Q = 2A - X VH J 9S ' - , 29 

((oo2-(qv) 2 )Vl FF ) 

valid for u), qv F <C A. 

From Eqs. (l20l - (!22j) . (|28|) and (|29|) one can obtain the complete BCS polarization tensor 

IP" (a;, q) in the vector channel. Due to conservation of the vector current the polarization 

tensor can be decomposed into the sum of longitudinal (with respect to q) and transverse 

components, where the longitudinal and transverse polarization functions are defined as 

U L = Iloo, n T = \ (dy - gig,/g 2 ) 1% with i,j = 1, 2, 3. 

8 



IV. FERMI-LIQUID EFFECTS 

We now turn to the Fermi-liquid effects. Our aim is to express the effective fields (J6]) via 
external fields. For this we can use the fact that the density current commutes with the 
bare interactions and does not need to be renormalized. The current connection with the 
external field A is given by the well-known relation 

5>,Cp,,) = £(j + M A ). (30) 

Inserting this into Eqs. fl3]) and (jlj) we obtain the time component of the effective field ((6]) 

in the form 

1 

CyQ 

Using also the continuity equation, qJL = up, we obtain the space components of the effective 
field: 



A o f = A o + ^rz F ^ P , (31) 



MM* u 



A e L " = A L -F 1 ^ r ^-5p (32) 

QPf c v 1 

A^ = A T - Fl ^j T (33) 

QPf c v 

Insertion of the effective fields @T) and (|33D into 5p = U L (A e ff - (u/q) A e / f ) and 

]t = —UtAj! after a little algebra results in the following: 

5 " = i-^ + ^/a + ^/sDepVn, ( A ° - ~ q AL ) ' (34) 



JT = ~- 7T7T, , r, /o^.-2 -2 _ 1TT A ^ • ( 35 ) 



Ik 

l-F 1 /(l + F 1 /3)v?c^Q-m T 
In obtaining this we have used the Landau formula relating the bare mass of a particle with a 

renormalized mass of a quasiparticle in the translation-invariant system M* jM = 1 + -Fi/3. 
The complete polarization functions LT^r relate the density perturbation and the density 
current with the external field as 5p = Ul (A) — {oo/q) Al), and jy = — UtAt, respectively. 
Then from Eqs. (|34|) and (135|) we obtain the complete polarization functions: 

UL = l-(F + S iF 1 /(l + F 1 /3))c v 2 Q-m L ' (36) 

Ut = i-F 1 /(i±i' 1 /3)v^ 2 Q-m T ' (37) 



in agreement with the results of Gusakov 18]. 



Equations (136]) and (157)1 completely describe the non-equilibrium behavior of superfluid 
Fermi liquids in the vector-linear-response region for not too high uj and q. Explicit expres- 
sions can be written using the following notation: 

, . , 1 1 / r. (cos 2 fl.,„ - g e'/El dn \ ,„„. 

7 (s ' h) = 2 U * s 2 -(«rfM^ d£/ ' (39) 

,-r v,l ' ! / fj co3 2 » qn (cos 2 qn - S 2 )f 2 /flg rf„ \ 

C (S ' h) = 6 + 2 U * s 2 -(cos 2 0,,„) £ 2 / g g 5^/ ' (40) 

" <S ' h) = U ^-(cos^ 2 /*?^/' <41) 

* (s - h) = (r V-Sifts^) • <42) 



where s — to/ (qvp), 



dn 1 , o -E D ,,„ s 

— cosh" 2 -E , (43) 



d£ p 2T 2T 



and we assume that the angle qn between the quasiparticle momentum p = pn and the 
momentum transfer q is fixed by the relation 

cos qn = h x sin 9 cos ip + h y sin 9 simp + h z cos 9 , (44) 

where the unit vector h = q/g = (h x , h y , h z ) defines the direction of the momentum transfer, 
and the unit vector n = (sin# cos<y?, sin# simp, cos#) defines the polar angles (9, tp) on the 
Fermi surface. 

Then the longitudinal polarization function ( 1361) can be written as 



Q(s,h) 



where 



UL ~ ^l-iFo + s^/il + F^QisM ' (45) 



Q (s, h) = n (s, h) + X (s, h) + 2a (s, h) ^ ^^P ^ • (46) 



V. SOUND-LIKE EXCITATIONS 

The pole of the density fluctuation propagator (fl"5|) at 

(F + S 2 Fx/ (1 + Fi/3)) Q (s , h; T) = 1 , (47) 

10 



defines the dispersion, s = so, of the " collisionless" collective mode, with u, qvF <C A. 
Equation ( HTj) with Q (s, h) as given in Eq. f H6|) generalizes previous results of Refs. (l9[ 



and 



20 



22J to the case of pairing caused by spin-orbit and tensor interactions. Therefore 
before proceeding to the detailed analysis of the sound propagation in the 3 P2 — 3 F2 neutron 
superfluid, we examine the obtained equations for the particular cases of the triplet-spin 
condensate in superfluid 3 He. 

Consider first the case of isotropic pairing. If the energy gap is isotropic, the angle 
integrals in Eqs. ( 13 8 1) -([421) can be performed by assuming the polar axis along the trans- 
ferred momentum. We then obtain 7 (s, h) = x (s, h) = and the longitudinal polarization 
function reduces to Eq. ( 1451) with 

^ / n / n 2a (s) C (s) , ,„ N 

q( '> = "M + *,(.)-<(.) ■ (48) 

I — I 

as obtained by Leggett |19j for the case of isotropic s-wave pairing. 

Consider now the case of anisotropic pairing at T = 0. In this case the quantities 
a (s, h) ~ 1/2 and £ (s, h) = 1/6 are independent of s and again isotropic, while 7 (s, h) = 
1] (s, h) = x ( s ) h) = 0. Thus for T = we find 

fl '= (itF,)..-. 8 -' S = 5 (1 + F ° )(1 + fl/3) - < 49 ' 



in agreement with the result obtained by Wolfe 21[ for an anisotropic phase of superfluid 
3 if e at zero temperature. Notice that the same was obtained also by Leggett [19] for the case 
of isotropic s-wave pairing. The pole at s = so corresponds to the first sound ("Bogolyubov- 
Anderson" mode) undamped at zero temperature. 

We turn now to the case of anisotropic p-wave pairing in liquid 3 He. The experimentally 
observable sound velocity in liquid 3 He is large, therefore it is traditional to calculate the 
sound dispersion in the limit s ^> 1. Expanding Eq. ( 146|) in powers of 1/s, one can obtain 
up to accuracy s -4 , 

«(^)4^ + K- <A(n)C084 ^ >+i ^5^)' (50) 

where 

/oo a 2 771 

de— tanh— . (51) 

11 



In obtaining Eq. ( )50|) we used the identity 

, e 2 dn f°° A 2 , F p 

de——— = l- de— -tanh-£ . (52) 

EldE p J El 2T y } 

Inserting Eq. fl50|) into Eq. (T4"T|) one can obtain the dispersion law for high-frequency sound 
in superfluid 3 He. Assuming F ~ Sq ^ 1 we find 



+ (3 + FO ( - <A (n) cos 4 tf nq ) + (MnW 2 g nq ) 2 ] , (53) 



in agreement with the expression derived by Wolfe 21]. 

We focus now on the sound propagation in the 3 P 2 — 3 F 2 superfluid neutron liquid which 
is expected to exist in neutron stars at supernuclear densities. First one has to specify the 
order parameter ([7]) for the particular case of neutron pairing. It is conventional to represent 
the triplet order parameter of the system as a superposition of standard spin-angle functions 
<&ji m , of the total angular momentum (j = 2,m,j) with partial amplitudes A^ m : 

D = ^A, m ^ m .(n) . (54) 

Irrij 

In our calculations we use vector notation which involves a set of mutually orthogonal 
complex vectors hi mj (n) defined as 

b, m . (n) = - (1/2) Tr (g&$ jlmj ) (55) 

and normalized by the condition (b*/ m '.b/ mj . > = $ii>S mjm '.- We will focus on the 3 P2 — 3 F 2 
condensation into the state with mj = which is conventionally considered as the preferable 
one in the bulk matter of neutron stars. In this case one has 

b 2 (n) = i (1 + 3n 2 ) 8 X + 2 (5n 4 - 2nj + l) S 3 , (56) 

where 5± = A 2 /A 2 and 5$ = A|/A 2 are partial contributions of the 3 P 2 and 3 F2 states, 
respectively, Si + S3 — 1. 

For < T < T c the behavior of il^ in the intermediate region of s depends essentially 
on the temperature. According to Eqs. (I3"8j) - (j4"2l) the imaginary part of the functions arises 
from the pole of the integrand at s 2 = (cos 2 # qn ) e 2 /E 2 This is Cherenkov's condition 

12 



which can be satisfied only if s < 1. Neglecting the narrow temperature domain where 
the imaginary part of polarization is exponentially small 19|, one can conclude that the 
well-defined (undamped) waves correspond to s > 1. Further we consider only undamped 
sound-like oscillations with s > 1. 
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FIG. 2. Zero-sound velocity as a function of temperature below the critical point, T < T c . 
The dimensionless velocity of zero sound, sq = u)o/(qvp) is shown vs the temperature parame- 
ter A (T) jT for different sets of Landau parameters. Case (a): Fq = 0.2, 0.6, 1; F\ = 0.5. Case (b): 
Fq = 0.2,0.6, lji^i = I. We present the curves for the wave propagation along the symmetry axis 
(0 q = 0) and in the perpendicular direction (0„ = 7r/2). 



The density-dependent Landau parameters entering the dispersion equation are not re- 
liably known. Therefore in Fig. [2] we present solutions to Eq. (|47p for several sets of the 
Landau parameters. The curves show the zero-sound velocity as a function of temperature 
parameter A (T) jT for a pure 3 P2 pairing (solid curves) and for the case of pairing into the 
mixed z Pi — 3 F 2 state (dashed curves). The plots are made for the sound propagating along 



the axis of the wave function of the condensate 



0) and in the orthogonal direction 



(0 q = vr/2). 

For the case of mixed pairing we have chosen A3 = 0.17Ai, in agreement with that found 
in realistic calculations by different authors (see, e.g., Ref. |6|). As one can see, a small 
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admixture of the 3 F2 state does not modify markedly the dispersion curves obtained for the 
pure p-wave superfluid. 

The sound waves are anisotropic. At fixed temperature the velocity of the sound grows 
along with deviation of the wave vector from the axis of the order parameter. The sound 
speed is maximal for orthogonal propagation. 

As regards the temperature dependence, immediately below the critical temperature T c 
the velocity of zero sound goes down, and the undamped wave disappears at some temper- 
ature T\ (0 q ) < T c when the sound velocity becomes smaller than the Fermi velocity, So < 1 
(although the mode with exponentially small damping can exist in some region below this 
temperature). Thus the undamped collective excitation may or may not exist at some tem- 
perature, depending on the values of F and F 1 . If it does, its velocity s (6* q , T) will always 
be greater than unity. 

VI. SUMMARY AND CONCLUSION 

Let us summarize our results. We have studied the linear response of a superfluid neutron 
liquid to an external vector field in the limit u,qvF <C A. The calculation is made for 
the case of 3 P2 — 3 -^2 condensate which is expected to exist in the superdense core of 
neutron stars due to spin-orbit and tensor pairing interactions. By analyzing the poles of 
the longitudinal response we have found the low-energy spectrum of sound-like collecive 
excitations caused by density fluctuations in the condensate. Previously the sound-like 
excitations were investigated for a triplet condensate caused by central pairing forces in 3 H e 
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22j . Our dispersion equation (|47|) represents a generalization of the above results to 
the case of pairing caused by noncentral spin-orbit and tensor interactions and naturally 
recovers previous results obtained for the case of central forces. 

The sound-like spectrum of a Fermi liquid substantially depends on the residual particle- 
hole interactions which are conventionally described by a set of Landau parameters. We 
have limited our consideration to the first two terms of this expansion. This approach can 
be considered as a model of the Fermi-liquid interactions, although there are indications 
that the higher- order Landau interactions do not affect the longitudinal response functions 

nn 

|26|,|27|. 

Unfortunately, the density- dependent Landau parameters entering the dispersion equa- 
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tion are not reliably known for an asymmetric nucleon matter (although, in principle, these 



can be evaluated theoretically 29H32J). Therefore we have studied solutions to Eq. fj47|) for 
several sets of the Landau parameters. We found that the sound waves are anisotropic and 
a small admixture of the 3 F 2 state does not modify markedly the dispersion curves obtained 
for the pure p-wave superffuid neutrons. At fixed temperature the velocity of the sound 
grows along with deviation of the wave vector from the axis of the order parameter. The 
sound speed is maximal for orthogonal propagation. 

Immediately below the critical temperature T c the velocity of zero sound decreases when 
the temperature goes down, and the undamped wave disappears at some temperature 
T\ (# q ) < T c when the sound velocity becomes smaller than Fermi velocity, (although the 
mode with exponentially small damping can exist in some region below this temperature). 
Thus the undamped collective excitation may or may not exist at some temperature, de- 
pending on the values of F and -F\. If it does, its velocity will always be greater than Fermi 
velocity. 

In our analysis, we have assumed that the axis of the order parameter is equally oriented 
everywhere. It is necessary to notice, however, that texture effects can orient different parts 
of the sample differently, and therefore give a range of frequency shifts which together appear 
as a broad line. The texture effects can be minimized by an external magnetic field, although 
the magnetic field has no consequence on the dispersion of the sound wave except that it 
serves to fix the relative orientation of the spin-orbital wave function associated with the 
order parameter, if the dipole interaction is taken into account. 

On the other hand the texture effects could play an important role in neutrino cool- 
ing of neutron star at the latest stage. Indeed, the sound wave can emit a neutrino pair 
through neutral weak currents while crossing the border where the axis changes its direc- 
tion. Neutrino radiation is possible also due to collisions of sound waves [131 ] . As already 
mentioned in the Introduction, the sound waves are known to play an important role also 
in a superffuid heat conduction when the transverse electron motion is strongly suppressed 
by a magnetic field [12J. Various applications of the results obtained in this paper will be 
considered elsewhere. 

All the results of this paper depend on the Fermi-liquid functions Fq, F\ which parametrize 
the normal Fermi liquid. The only way to evaluate these functions for a superdense asym- 
metric nuclear matter is to estimate their values from first-principles calculations. Although, 
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in practice, such work is in progress |29l432| . the complete information on the Landau param- 
eters for a neutron matter is not available at the moment. Only the well-known conditions 
[171 ] of the matter stability with respect to long-wave perturbations can be used in order to 
limit the Landau parameters. 
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